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ABSTRACT 

We present a non-parametric measure of the ellipticity and the alignment of the 
dark matter halos in Abell 901/902 supercluster. This super-cluster is a system of 
four separate peaks in a 0.5° x 0.5° field of view. We map the mass distribution of each 
individual peak using an improved version of Particle Based Lensing (PBL) and measure 
the ellipticity of the dark matter halos associated with two of the peaks directly from 
the mass map and by fitting them to a singular isothermal ellipse. The parametric and 
non-parametric measurements are consistent for A901b while the position angle for the 
Southwest Group is different for the two techniques. We account for this discrepancy 
to substructure present in the Southwest Peak. We estimate an axis ratio of 0.37 ± 0.1 
for A901b and 0.54±{j;gjj for the Southwest Group. 

Subject headings: methods: statistical, methods: analytical, methods: data analysis, 
gravitational lensing , galaxies : cluster s 



1. Introduction 



Clusters of galaxies are the largest virialised struc tures in the universe. Their mass distri- 
bution reveals propertie s of the primordial density field (jSheth &; Tormenlll999l : iFrenk et aJj[l990j ; 



Press fc Schechterl 119741), and their mass function is sensitiv e to predictions of ACDM cosmology 



([Francis et al,l |2009| ; iBartelmann et al,l 120061 ; lEke et al,l Il996l ) . As suggested by their low gas frac- 
tion, the baryons in clusters have little influence on the mass distribution beyond the core. Thus, 
clusters can be well approximated as dark matter halos. These halos produce spectacular arcs and 
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mult iple images (jZitrin et alj|2009l ; IZitrin fc Broadhurstll2009l ; iLimousin et al.ll2008l ; iGoldberg et al 
20091 ) and cause weak distortions of background galaxies. 



Weak lensing inversion techniques study this distortion pattern to infer properties of the 
mass distribution of th e cluster independent of its physical/dynamical state (IClowe et al.l 120061 : 



Okabe fc Umetsu 



200 



Weak lensing has been used to measure dark matter density profiles 



in clusters (|Broadhurst et al.1 l2008l . 120051 ; lUmetsu fc Broadhurstl 120081 ; iMandelbaum et al 



2008) 



to test predictions o f N-body simulations of the standard ACDM model (jSpergel et al.l 120071 ; 
Komatsu et al.l l2009i ) . However, weak lensing has systematic and statistical errors. We have to 



understand and reduce these errors in order constrain the mass, size and shape of galaxy clusters. 



Gravitational Lensing mass reconstruction of clusters have been studied for two decades ( 


Okabe et al. 


200fl|; 
2000|; 


Kine & Corless 


20071: 


Pedersen & Dahld 20071: Ireens et a 


J2002I: 


CvDriano et alJboOl: A 


/an Waerbeke et al. 


Williams et al. 


19991: 


AlleJ 


1994 


Wu & Hammer! Il995l: r . 


?vson et al.l 199ol: Oeuri et al 


20091: 


Bardeau et alJ2007:] 


Bradac et al. 


2006: 


Broadhurst et al.ll2005: 


Clowe & Schneider 


2001: Bacon et al. 


2006) and many different methods have been developed in the process ( 


Kaiserlll995l: 


Seitz & Schneider 


1995: 


Bartelmannlll995; Schneider & Bartelmann 1997; Seitz & SchneiderlboOl). Almost all of them 



perform reasonably well in detecting massive peaks (primarily associated with the Brightest Clus- 
ter Galaxy) in the field of view. This is because massive peaks produce significant lensing signal 



of observations has lead to extensive research on lens mass reconstruction techniques ( 


Seitz 


et a . 


1998 




Geieer & Schneider 


1999: Marshall et al. 


2002; 


Lombardi & Bertin 


199 


J; 


Dieeo et al. 


2007: 


Bradac et al. 


20051: Merten et al. 20091: 


Bradac et al.l 


2009; 


Mandelbaum et al. 


20091). 



In this paper we have improved Particle Based Lensing (PBL) (|Deb et al.ll2008f ) by smoothing 
the ellipticities prior to mass reconstruction and evaluated the covariance of the resulting mass 
distribution. PBL is a mass reconstruction technique where the mass (or potential) is calculated 
at the location of each image. There is a weight function associated with each image with a kernel 
having a width that varies according to the local number density. 

We aim to measure the ellipticity of the dark matter halos of the super cluster Abell90 1/902. 
Abell 901/902 is a complicated system with several sub-clusters spread around a field of view of 
0.5° x 0.5°(~ 5x5 Mpc 2 ) at a redshift of z = 0.165. This field was surveye d to study how 
galaxy evolution gets affected by the density of the environment (jGray et al.ll2009i ). The mass and 
morphology of the dark matter halos of these galaxies play an important role in their response to 
surrounding environment. Since this field of view has three clusters and a group of galaxies it is a 
very good laboratory to study large scal e struc ture and its influence on galaxy transformation. The 
weak lensing analysis (jHeymans et al.l (|2008l ). hereafter H08) of this field reconstructs the large 
scale dark matter for this super cluster. The mass estimates of each peak were made using NFW 
profiles and non-parametric maximum likelihood methods. Comparison with the light profile has 
shown that the substructure in the dark matter peaks are closely followed by the substructure in 
galaxy groups. 
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The paper is organized in the following way. In § [21 EHU we discuss the technical aspects of mass 
reconstruction and noise covariance analysis using PBL. In § [5] we describe measuring ellipticity 
of the dark matter halo using parametric and non-parametric techniques. In § [6] we give a brief 
description of the data and in § [7] we report the estimated ellipticity for the individual peaks of 
the super cluster Abell901/902. In the last section we discuss our results and possible directions of 
future work. 



2. Weak Lensing Reconstructions 

Weak lensing is a statistical measure of small distortions in background galaxies caused by a 
cluster. The distribution of intrinsic elli pticities in so urce galaxies follows a Gaussian distribution 



with zero mean and a width of 0.2 ~ 0.3 (|Mellierlll999l ). In order to extract a lensing signal from the 



weak distortion of image galaxies several of them have to be averaged to smooth out the intrinsic 
noise. The ellipticities of the background image galaxies are the weak lensing observables. In the 
linear regime the ellipticity can be approximated by the shear of the lensing field. The shear has 
two components given by, 

7 = 7 1 + i^ 2 (1) 

The shear 7* and the dimensionless surface mass density k are linear combinations of the second 
derivatives of the potential in angular coordinates, 

^1 _ 1p,xx — 1p,yy ^ 

7 2 = ^,xy (3) 

(4) 



Y,xx V" ',yy 



2 

where ip^xx^.xy, ■•• refers to jp^:, dxdy i •■■ ^ n ^ ne weak lensing regime the ellipticity induced by the 
lens is given by, 

e '= 9 ' = (T^) (5) 

and in the strong lensing regime the induced ellipticity is given by the inverse conjugate of the 
reduced shear. In this paper we will be addressing the weak lensing regime only. The above 
expressions are written for a source plane at infinity. The redshift dependence for convergence and 
shear are included by, 



k(0, z) = Z{z)k{6) Y(9, z) = Z(z)y(9) (6) 
where Z(z) is the redshift weight function for background images given by, 

Z(z) = ^ (7) 
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2.1. Non-parametric mass Reconstruction Techniques 



Non parametric weak lensing mass reconstruct ion techniques are broa dly categorized into two 
categ ories. There are direct t echniques like KS93 (jKaiser h Squireslll993l ) and Finite Field Inver- 



sion (Seitz &: Schneider 



2001 



j) where a linear relation is assumed between the measured ellipticity 
and the convergence. The advantage of these methods is that they are very fast and useful for 
testing the data. The disadvantage is that the linear techniques reconstruct the non-linear regimes 
less accurately. Hence the inverse techniques become important as we approach the strong lensing 
regime. In these methods a likelihood function is written between the observed ellipticties and the 
reduced shear and maximized iteratively to obtain the best solution. 



C = exp 



E 



. e 

mn \ m 



c 



-1 



1-KmOm) J ^ mn \ £n 1-Kn(f„) 



(8) 



where e m \ are the observables. Using finite differencing or Particle Based Lensing (PBL) the 
derivatives of the potential can be written as, 

^ = D^ m (9) 

where represents derivatives of the potential, v is the order of the derivative. For first deriva- 
tives v = x, y, for second derivatives v = (xx, xy and yy) and so on. The likelihood function given 
by Equation [5] is linearized and written as a function of Dn m , the data correlation, the potential 
and the constraints at each step of the minimization and the potential is solved iteratively at the 
maxima of the likelihood. 



3. Method optimization 

Weak lensing is a statistical measure of small distortions in background galaxies caused by a 
cluster. In order to extract a lensing signal from the weak distortion of image galaxies several of 
them have to be averaged to smooth out the intrinsic noise. Smoothing is an integral part of any 
problem with noisy data and low signal-to-noise ratio. The error budget of weak lensing is controlled 
by the scale at which the data is smoothed. Optimizing this scale is necessary independent of the 
technique that is used. In this section we will introduce smoothing, describe the inversion technique 
used to create mass maps from measured ellipticities and lay down a foundation for the choice of the 
optimal smoothing scale that will produce minimum reconstruction errors in the recovered mass. 

In order to give equal weight to all ellipticity measurements around a single image we use an 
azimuthaly symmetric smoothing function. We choose a normalized gaussian for this purpose. 

where e l represents smoothed ellipticity field and Q is the smoothing function. Here m, n represent 
the background image index. Smoothing of data prior to mass reconstruction has been done by 
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several groups (|Bartelmannlll995l : ISeitz & Schneider! Il996l . Il998l . l200lh . the covariance due to this 
smoothing is given by, 

Cmn — QkmQkn&n (H) 

Here a n is the noise due to intrinsic ellipticities. This covariance between the constraints is used 
in the likelihood analysis in Equation 



3.1. Fitting the error 

We have defined a smoothing function and described the inversion procedure to create mass 
maps. The input parameter to this method is the smoothing scale. In this section we fit the 
weak lensing reconstruction error as a function of the smoothing scale, measurement error, number 
density and the length scale at which structure can be resolved. 



3.1.1. Toy Problem: A Sine Wave 

In order to understand how errors propagate into the reconstructed mass we will study a simple 
problem. The convergence k for this field is defined by, 

K = Asmi— J (12) 

where A is the wavelength, it represents the scale at which structure can be resolved. We do a 
weak lensing reconstruction of this wave and determine the scale at which the reconstruction error 
is minimal. We generate the mock data by adding noise to the shear (given by Equation I12j) drawn 
from a Gaussian of width <r e . This noise is varied over a wide range to determine a fitting formula 
for the error in the reconstructed k. We define the error that we are trying to measure as the 
difference between the model (Equation. [T2|) and the reconstructed mass at every image location 
(for particle based approaches like PBL) averaged over the entire field of view. This error is given 
by, 

K rue ^reconstructed Y) (13) 

Here it is important to remember that while reconstructing k from data we will not be able to 
use Equation [13] as our goodness of fit since we do not know the underlying distribution of matter. 
In that case, the minimum % 2 will determine the best reconstructed mass. 

We smooth the ellipticity field using a normalised gaus sian of width C, The c o variance in- 



trodu ced by smoothing data has been thoroughly studied by lLombardi fe Schneider! ([20031 . 120021 . 



2001 



). We will construct an empirical form that fits Equation ([ 13f) by considering the effects of 
smoothing on mass reconstruction. The aim of smoothing is to average out the noise in the data, 
however, in this process we also wash out any structure smaller than the smoothing scale. We take 
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into account these two effects and fit them to the error in the reconstructed convergence from the 
mock data. 

W)~ °\V (CM) 2 +Co (14) 

Here n is the number density per unit wavelength of background galaxy images. The first term 
represents the second order bias due to smoothing of small scale structure. The error introduced by 
this term is proportional to the square of the smoothing length. The second term is the contribution 
from external error, in case of weak lensing this is the error due to the intrinsic ellipticities of the 
background galaxies. <t 2 /(k 2 ) represents the noise-to-signal ratio. The contribution from this term 

is inversely proportional to n, implying that the external errors get washed out as the number 
density of images increases or as the area of smoothing increases. 

Figure [1] shows a fit for equation [TH to the observed error for increasing amount cr 2 . We have 
plotted two cases, the solid (dashed) line represents the fitted value and the triangles (squares) 

2 

represent error for ^§y = 1(5) measured from the reconstruction. Each data point (triangles 
and squares) on the graph represents a different mass reconstruction corresponding to a different 
smoothing scale and noise. As is evident from the plot Equation Q3] is valid around the minima of the 
curves which defines the optimal smoothing scale for mass reconstruction. This happens because 
Equation [T5] was constructed by considering the first term in the Taylor expansion with non-zero 
contribution to errors. As we move away from the minima higher order terms start dominating. 
Since we are interested in computing the smoothing scale at which the error is minimum we do 
not need the higher order terms. This is a demonstration with a toy model, we generalize this 
form and use Equation [T5] to determine the scale at which we smooth ellipticities prior to the x 2 
minimization. 



Covariance 



In this section we will compute the covariance due to smoothing in the reconstructed mass, 
which will be used to compute the ellipticity of the mass distribution. T he covariance in the likeli- 
hood function has been studied extensively in the context of cosmic shear. lEifler et al.l (]2008bl ) have 
studied the covariance between cosmological p arameters to determ ine its effect on the likelihood 
analysis for cosmological parameter estimation. lEifler et al.l (|2008al ) uses the covariance among the 
data when comparing the information content in aperture mass and two-point correlation func- 
tion. However, not much importance has been given to the r ole of the covariance in cluster cluster 
mass reconstructions with three exceptions of iBridle et al.l (|1998l ) who derive an analytic expres- 
sion for the cqyarian ce and compare it with the covariances from Monte-Carlo simulations. Also 
van Waerbekd (|2000l ) have shown that a maximum like l ihood lensing mass reconstruction has a 
noise that follows a gaussian random field. iMerten et al.l (|2009i ) have also used covariance between 
the mass bins for computing mass maps. 



2 Z " 4 

a O a e /(™C ) + a lC + a 2 
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Fig. 1. — Error vs the smoothing scale for two levels of noise. The triangles and the squares 
represent the errors measured from mass reconstruction done for different values of noise and 

2 

smoothing scale. The solid line and the dashed line represent the fitted values for = 1,5 given 

by Equation \TQ here £ = 4 . The minima in this plot represents the ideal smoothing scale for 
the given noise. As expected decreasing the signal-to-noise for the mass reconstruction shifts the 
minimum towards higher smoothing scale. 

We have defined the covariance due to smoothing of external error in § [3] in the ellipticities. 
The covariance of any linear function f{9) sampled at positions 9 n between any two positions 9 n 
and 6 m is given by, 



For the sake of simplicity of notation in the rest of the paper we will denote Cov(/; n , 9 m ) = Cov4 m , 
where n, m represent either the total number of image galaxies in case of PBL or the total number 
of grid cells in case of a grid based method. We will derive the expression for the covariance 
matrix in the linear regime. This calculation is applicable to any technique where the derivatives 
of the potential can be expressed as a matrix times a potential. Hence it is applicable to PBL and 
finite-differencing techniques. In the very weak lensing regime where k<1 the ellipticty can be 
approximated by the shear, 7, 



where i=l, 2 represents the two components of ellipticity and shear. In this case we can write a 
likelihood between the observed ellipticty and the measured shear. Using Equation we can write 
the likelihood as, 




(15) 



(4) = < 



(16) 
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C = exp 



(17) 



where G l is matrix that relates the potential to j 1 , C is the covariance due to smoothing and e l is 
the smoothed ellipticity field. Using Equation 

ryxx I jjyy 

K = (18) 

jjxx _ jjyy 

G 1 = - (19) 

G 2 = D xy (20) 

Setting the first derivative of the log of the likelihood to zero we get, 

^(G^C^GV = ^(GYC- 1 ? (21) 

i i 

From this equation we can solve for the potential, once we know the potential there is a linear 
relationship between the measured ellipticities and k. In this case the estimated k is given by, 

k = km~ x ^y&fc- 1 ? = v*e (22) 

i i 

where K is matrix that relates the potential to k, M = ^2 i (G l ) T C~ l G l and 

V 1 = KM~ 1 (G i ) T C~ 1 (23) 

The covariance in k follows Equation [15j This can be re- written as, 

Cov K = FCov e V T (24) 

where Cov e is the covariance in the ellipticity. There are several effects that contribute to this 
covariance. As discussed in § 13.11 one term is due to smoothing of small scale structure and the 
other is due to averaging the error. We have derived the effect of smoothing on errors due to intrinsic 
alignment of galaxies and used it for the maximum likelihood analysis. Here we are going to derive 
a more general expression for the covariance in the final reconstructed mass, The covariance in 
measured ellipticity is given by, 

Cov mn = (( e m ~ 7m) {e n - %)} = 5 mn al (25) 

where 7 P is the true shear for the lens. After smoothing this becomes, 

Cov mn = ((£m ~ im)(£n ~ fn)) (26) 
— QmpQnpip^p ~\~ Tm'Tn) ^n^m "fmQnp'fp '~fnQmp'Yp 

Let us define, 

T?n — QmpTp (27) 
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Using this notation we can write the above equations as, 



Cov mn = QrnpQnp<?l + ((7m ~ 7m) {in ~ In)) (28) 

The first term is same as Equation QTJ Since we do not know the true shear we replace 7 with the 
reconstructed shear. Replacing this in Equation 1241 we have an expression for the covariance in the 
linear regime. These relations are summarized in Table [TJ 



M 


X 


y = mx 


Cov y 


G 1 




71 


GWCov^GW T 


G 2 




72 


G^Cov^'G^ T 


K 


Y> 


K 


KCov^K T 


Q 


7 


1 




Q 


£ 


£ 


Qaa T Q T + {(i~ m - lm){ln ~ In)) 


V 


£ 


K 


FCov £ V T 



Table 1: A summary of the relation between various matrices defined in§|H The last column gives 
the covariance among the observable in the third column. The last row is the final expression for 
the covariance in the reconstructed K. 



5. Where is the information in Cluster Lensing? Cluster Ellipticity 



It has long been establ ished from s i mulat ions (jjing &: Sutol [2002j; iRahman et al.l 120061 ) that 
galaxy clusters are triaxial. lOguri et al.l (|2005l ) have used triaxial dark mat ter halos to study the 
steep mass profile of A1689 and arc statistics using semi-analytic models (jOguri et al,ll2003l ) for 
triaxial dar k matter halos have ex plained the abundance of gravitationally lensed arcs for a sample 
of clusters. ICvpriano et al.l (|2004h have measured the ellipticity and the position angle of a sample 
of X-ray selected clusters parametrically using gravitational lensing and found that the dark matter 
halos were aligned with the brightest cluster galaxies. In this paper we will compute the ellipticities 
of dark matter halos non-parametrically. This model-free estimation is done by calculating the 
quadrupole moments of the mass map using the noise matrix derived in the previous section. 



There is a strong dependence of ellipticity on amplitude of mass fluctuations as (|Ho et al 



20061 ). A higher value of a$ indicates that clusters formation has started earlier and hence the 
measured ellipticity of clusters in the local universe would be lower. Clusters are formed through 
hierarchical merging of smaller dark matter halos. Thus at their infancy they have more infalling 
matter and are more elliptical. As they virialise they become more and more spherical. Thus we 
expect clusters at higher redshift to be more elliptical than low redshift clus ters. This has been 
confirmed by measu ring higher order moments of the X-ray gas distribution (jjeltema et al J 120051 : 
Buote fc Tsailll995l ). Following the procedure that we lay down in this paper we can measure the 
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ellipticity of cluster halos from lensing for a large sample of clusters distributed in redshift. This will 
make the contribution of the errors due to projection smaller. Since gravitational lensing probes 
the projected mass of the lens it is very difficult to constrain whether the halo is oblate or prolate. 

We measure the shape of the lens by calculating the second order moments of the mass distri- 
bution, this will give us the eccentricity and the position angle for an elliptical mass distribution. 
We will also fit the ellipticity of the lens parametrically, for a truly elliptical lens this will give a 
very good description of the shape of the underlying dark matter halo. 



In the previous sections we have derived the covariance of a mass map in the linear regime. The 
super cluster Abell 901/902 is a sub-critical cluster, hence for the error analysis we assume linearity. 
It is important to note that linearity is not assumed for doing the mass reconstruction. Since we 
have a correlated mass map with a covariance that is singular, we will perform a singular value 
decomposition as explained in § |A] and consider a few eigen modes that are significant. We then 
transform to a basis where the eigenmodes are independent. This is done by the transformation, 



where U is the orthogonal matrix from singular value decomposition defined in Equation I All 
The k"s are independent. We will express all quantities in terms of n'. Here (k) is the mean density 
of the field, it is calculated as follows 



The ellipticity of the dark matter is defined in terms of the quadrupole moments. The simple 
definition for quadrupole moments for the k field is given by, 



We know that k is correlated and hence we write the above expression in terms of n 1 and 
inverse weight it by s, where s is defined in Equation IA11 The errors in the Quadrupole moments 
are given by the Octopole moments and are calculated in a similar fashion. The ellipticities of the 
lens are defined by, 



5.1. Measuring Cluster Shapes Non-parametrically 




(29) 




(30) 




(31) 




(32) 
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C J 2 



Ixx + lyy + {Ixxlyy ^xy) 2 



(33) 



5.2. Measuring shapes parametrically 

We will fit singular isothermal ellipse to the dark matter peak s to constrain the sha pe of the 
dark matter halos. The convergence and the shear for this profile (jKormann et al.lll994l ) is given 
by, 



K = 7 



19e 
2 6 



f [cos((f> — a) + f 2 sm((f) — a)] 



(34) 



where Oe is the Einstein radius given by, 



9 E = 4vr^^ 



c 2 D„ 



(35) 



and a v is the velocity dispersion of the cluster. / is the axis ratio and a is the position angle. This 
relationship is degenerate with a = a + tt/2 and / = 1//. The complex ellipticty is related to the 
to axis ratio via, 



1-/ 
l + f 



(36) 



This model is fitted by minimizing a x 2 of the form 



X 



E 



(37) 



where i = 1, 2 the two components of ellipticity, e l m is the ellipticity of the mth galaxy and g l m is 
the reduced shear an d o> is the error in the ta ngential ellipticity. This form of the mass model has 
been investigated by King &: Schneider! ()200lh and applied to a sample of X-ray luminous clusters 



(|Cvpriano et al.l 12004 ) 



6. Data 



Th e shear catalog f or thi s cluster is generated following algorithms described in lKaiser et al 



(1995): 


Hevmans et al. 


Hevmans et al. 


( 


2005). 



modelling of the PSF, the stars were chosen to maximize the signal-to-noise of the stellar ellipticity 
function and minimize the temporal variation. The charge transfer effieciency of the ACS has 
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degraded over the years, this is corrected using the methods proposed by iRhodes et al.l (|2007l ). 
Furthermore magnitude cuts are applied to the galaxy catalog to eliminate cluster members and 
foreground galaxies. The magnitudes are chosen such that 23 < mpQ^w < 27.5 and the signal to 
noise for each galaxy is chosen to be > 5 with a size greater than 3 pixels. The total sample size 
is ~ 60000 galaxies with 65 galaxies per square arcminute. For Abell901/902 the redshift for 90% 
of the source galaxies were not known. However, the redshift weight function does not have a very 
strong dependence on the redshift for z s > 1 and zi = 0.165. Hence we assume the sources to be at 
a single redshift of 1.4 which is the estimated median redshift of the background sources. A more 
detailed description of the data and tests for systematics can be found in H08. 



7. Results 

We have outlined a recipe for choosing the smoothing scale and using PBL for mass reconstruc- 
tion. We will apply this technique to the Abell 901/902 fie ld of view. The dark matter peaks in 



this field of view has been detected bv lHevmans et al.1 (|2008l ). We will zoom into each of the peaks 



and reconstruct a mass map and an error map for the peaks. For these peaks the signal-to-noise 
ratio ~ 0.3, using this value in Equation 1141 we get an approximate smoothing scale of 0.5'. We 
reconstruct the dark matter distribution of the four peaks. First we will study A901b, it is the most 
compact peak with a single halo. The minor to major axis ratio and position angle is reported in 
Table [2j The ellipticity of the peak is 0.45to'{o ano - the position angle is 90° suggesting that it is 
elongated vertically. The top left panel of Figure [2] is a map of A901b. 

The bottom left panel of Figure [2] is the dark matter reconstruction of the Southwest Group. 
This peak clearly shows that it is not spherically symmetric, it is elongated in the north west 
direction with significant substructure suggesting that the group is not completely virialized. We 
measure an axis ratio of 0.3 ± 0.07 and a position angle of 120° ± 4.8°. Since the dark matter 
distribution has multiple maxima an elliptical dark matter halo may not be the best description 
for the Southwest Group. A901a has two distinct peaks. This is clear from the dark matter map 
of the top left panel of Figure It has non-zero quadrupole moments, however, it is not possible 
to describe it using an ellipse. 

If we compare Figure 12131 with H08 we see that the dark matter distribution is very similar. 
However, in this analysis we have smoothed the data prior to reconstruction and included the 
covariance due to smoothing in the x 2 minimization and in calculation of covariance of the final 
mass reconstruction. This makes the errors of the reconstruction well understood. The error maps 
in the right hand panels of Figure [21 [3] is computed by taking the square root of the diagonal of 
the covariance matrix of k. The central peak of A901b is detected at 7a significance, the two 
sub-peaks of the Southwest Group is detected at 4cr. The central peak of A901a is a 5a detection 
the secondary peak is a 2a detection. The peak of A902 is detected at 4c. The mass measured 
within one arcminute of the center of each peak is listed in Table [21 We measure the ellipticity 
of A901b and the Southwest peak since these two peaks are relatively less disturbed. A901a has 
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Fig. 2. — Mass reconstruction of dark matter halos in A901/902. Top Left Panel: A901b Top 
Right Panel: Error Map of A901b. Bottom Left Panel: SouthWest Group. Bottom Right Panel: 
Error Map for the Southewest Group. There is an artificial increase of error towards the extemities 
of the maps where the convolution kernel steps over its hard edges. A901b is a compact dark matter 
halo, the peak is detected at 7a and the Southwest Group has significant substructure with two 
sub-peaks detected at 4<r . We measure the ellipticity of these two peaks. 
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two distinct peaks and A902 is very disturbed even in the outer regions. Hence it is difficult to 
represent these peaks with ellipses. 



7.1. Parametric Fitting 

From the non-parametric reconstruction it is clear that A901a and A902 are double peaked 
systems, while A901b and the SouthWest Peak are closer to a single halo with substructure. We 
have fitted this data to a singular isothermal ellipse described in § 15.21 We fit a single halo centered 
on the BCG in each cluster. Since the aim of this study is to measure the ellipticity of the lens 
we do this fit for the SouthWest Peak and A901b only because A901a and A902 have irregular 
structure and hence cannot be modeled as an ellipse. The field of view contains four distinct peaks 
we consider a patch of 178 square arcminutes around each peak. This ensures that the shear signal 
is not contaminated by the other dark matter halos. The constraints on the velocity dispersion, 
ellipticity and position ang l e is li sted in Table [2j The constraints on the velocity dispersion are 



consistent with iGrav et al.l (J2002J) within error bars. We estimate the ellipticity of A901b to be 
0.39 ± 0.09, and a position angle to be 90° implying that the dark matter halo is elliptical in the 
vertical direction. The dark matter map of the southwest peak is clearly not spherical. This peak is 
well fitted by the singular isothermal ellipse and the ellipticity of this peak is 0.4^q \q. In Figured] 
we have plotted the joint 1(2) — a error probability distribution between the axis ratio and the 
position angle. It is not possible to constrain the ellipticity of A901a and A902 with an elliptical 
mass model. In fact a parametric fit is consistent with the spherical model with high error bars. It 
is clear from the mass reconstruction that the peaks are disturbed. As a matter of fact A901a has 
two distinct peaks, the second peak coincides with an infalling X-ray group. A902 also has another 
galaxy group in the background at a redshift z = 0.46. 



7.2. Comparison between dark matter and light distribution 

We have measured the ellipticity of the light distribution by measuring quadrupole moments 
of the cluster member galaxies weighted by their stellar masses. The results are listed in Table EJ 
The light distribution and the dark matter distribution is coincident for A901b. In the Southwest 
peak the light distribution is coincident with one of the sub-peaks. The light distribution is less 
elliptical compared to the mass distribution. The results from the non-parametric reconstruction 
of the Southwest group indicates that the major axis of the light distribution and the dark matter 
distribution are not coincident. This is because the major axis of the Southwest peak is influenced by 
the smaller sub-peaks whereas the light distrbution has a single peak. Figure [5] shows a comparison 
of the light and dark matter distribution. The background map represents the stellar mass of the 
cluster member in units of Mq. The contours represent the dark matter distribution. 
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Fig. 3. — Mass reconstruction of dark matter halos in A901/902. Top Left Panel: A901a Top Right 
Panel: Error Map of A901a. Bottom Left Panel: A902. Bottom Right Panel: Error Map A902. 
A901a has two distinct peak, the central peak is detected at 5<r and the secondary peak is detected 
at 2d. A902 is reasonably disturbed with the central peak detected at 4c These two peaks are 
not representative of elliptical dark matter halos. 
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a a 

Fig. 4. — The two panels are the error ellipses for A901b and the Southwest Peak derived from 
the parametric modeling, described in § 15.21 . The y-axis / is the axis ratio and the x-axis a is the 
position angle defined in section § 15.21 in radians. For both plots we have plotted the la and 2a 
contours. Both peaks have non-zero ellipticity at 2a level. 




Fig. 5. — A comparison of the light distribution vs the dark matter distribution. The colors 
represent the light distribution and the over-layed contours represent the dark matter distribution. 
Left Panel: A901b. Right Panel: Southwest group. The units of the light distribution is Mq. 

7.3. Comparison between parametric and non-parametric results 

We have computed the ellipticity and the position angle for A901b and the Southwest Peak 
parametrically and non-parametrically. The results for A901b using both methods are consistent. 
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The position angle for the Southwest Peak is inconsistent between the measurements. This peak 
has an irregular mass distribution, hence an elliptical mass model does not provide a full description 
of its shape. The va l ues of the measured ellipt icity are very close to that expected from simula- 
tions (|Ho et al.l 120061 ; iBailin fc Steinmeta 120051 ) and previous observations of dark m atter galaxy 
halo ellipticities from the Red Sequence Cluster Survey(RCS) iHoekstra et al.l (|2004j ) survey. The 
measured ve locity dispersion for the peaks are consistent with the results from ground based in- 
vestigations (jGray et al.1 12002 ). The error in the minor to major axis ratio is quite high since weak 
shear is a noisy estimator of the dark matter distribution, however it is the only way to uniquely 
measure the dark matter halo shape. 



8. Discussion and Future Work 



In this paper, we have measured the projected ellipticity of dark matter halos non-parametrically. 
This is done using an improved PBL by including smoothing of the ellipticity field. We have also 
calculated the covariance of the resulting mass distribution making the errors of the reconstruction 
well understood. We have applied this technique to the super cluster A901/902 and reconstructed 
each of the peaks individually and measured the ellipticity of A901b and Southwest Peak from the 
PBL reconstruction and using parametric models. The other two peaks A901a and A902 have a 
lot of substructure and cannot not be modeled as ellipses. We have not considered the line-of-sight 
ellipticity of the d ark matter ha l os of Abell 901/902. This is because lensing probes the projected 
mass distribution. ICorless et al.l (|2009i ) have fitted triaxial NFW to galaxy clusters with high errors 
on the concentration and mass suggesting that more information form strong lensing, X-rays and 
Sunyaev Zeldovich effect is required to constrain p arameters pertaining to the line-of-sight. A joint 
analysis using X-rays and lensing has been done bv lMorandi et al.l (|2009i ) to constrain the ellipticity 
along the line of sight. 

In future we will measure the ellipticity o f the projected dark matter halos for a large sample 
of clusters and compare them to simulations ([Hopkins et al.ll2005l ). The error introduced due to 
the line-of-si ght ellipticity f or a l arge sample will be insignificant. A similar study has been done 
in X-rays bv lJeltema et al.l (|2005l ) by studying the ratio of the higher order moments to the zeroth 
order moments. 
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A. Inverting the Covariance matrix for % 2 minimization. 

In order to do a % 2 minimization the inverse of the covariance matrix is important. Before 
doing the inverse it is important to understand the properties of the covariance matrix. It is essen- 
tial to remember that every image location is not independent. This is because we are smoothing 
the data, this makes the number of image positions larger than the numb er of 'resolution units ' 



for a given smoothing scale. T his situation has been encountered before ()Pan fc Szapudil 120051 ; 



Gaztahaga &: Scoccimarrd 120051 ) when the data vector was larger than the number of indepen- 
dent information content. This makes input data degenerate, i.e for a smoothing scale consider- 
ably higher than the interparticle separation two positions spatially close to each other essentially 
represnt the same information. This makes the covariance matrix singular. The number of in- 
dependent components of the covariance matrix is inversely proportional to the area under the 
smoothing kernel. 

We invert this singular matrix by using the Singular Value Decomposition (SVD) implemen- 
tation of Tikhonov regularization. The traditional SVD matrix inversion is given by, 



M _1 = VSU T (Al) 

Here U and V are orthogonal matrices, since C is symmetric U=V. S is a diagonal matrix , the 
diagonal elements are given by {1/s}. Here s represents the eigenvalues of C. In case of a singular 
matrix some of the eigenvalues are zero. Usually s is written in descending order and the first I 
non-zero eigen values are used for the matrix inversion. For the n — I eigen values 1/s is replaced by 
zero. In real cases it is commonly seen that the eigen values are not zero, rather they are numerically 
very small, and hence dominated by round-off error. These problems ar termed as ill-conditioned 
problems. In order to deal with this situation a truncated SVD is used to do the matrix inversion, 
i.e eigenvalues below a certain threshold are considered to be zero. Choice of this threshold is 
dependent on the particular problem. If the singular values of a matrix can be distinguished from 
the non-singular ones in a well-defined fashion then the choice of the threshold becomes simple. 
This is the case when there is a substantial gap between the largest singular eigenvalue and the 
smallest non-singular eigenvalue. However, in this particular problem where the covariance is due 
to gaussian smoothing of the data there is no such distinction. As a matter of fact the eigenvalues 
smoothly asymptote toward zero. Hence the problem does not direct us toward any obvious choice 
of the threshold value. Since the problem is severly ill-posed we use Tikhonov regularization(ref) 
to invert the covariance matrix. This is equivalent to 



M" 1 = VfU T (A2) 

where f m = 2 S T s are the filter factors. For s m » a, f m = l/s m and for s m « a, f m = 0. 
The presence of the regularization ensures that there is a smooth transition between f m = l/s m 
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and f m = instead of an abrupt cutoff threshold. The regularization parameter a is given by, 



As a decreases more and more eigenvalues are included in matrix inversion. As the smoothing 
scale is increased the area becomes higher, a becomes higher and the number of modes used in the 
matrix inversion becomes lower. 



In the previous section we have evaluated the covariance for the reconstructed k in the linear 
regime. Inverse techniques like PBL are often able to reconstruct the semi-strong regime of weak 
lensing clusters within a couple of iterations. Hence we write down a general formalism in the 
semi-strong region. 

In the second iteration the potential is given by, 



the correct solution for the potential both components of the ellipticity will be summed, for the 
sake of simplicity of notation we have not included it in Equation IBll We assume that we are in 
the semi-strong region and k < 1, hence the modelled ellipticity id given by, 



a = C( 2 



(A3) 



B. Generalization to the Non-linear regime 




(Bl) 



where the modelled ellipticity from the first step of the minimization. In order to get 




(B2) 



The reconstructed k in the second iteration is 



g 



;iven by, 



(B3) 



The expectation value for the ellipticity is given by, 




(B4) 



The covariance in k after the first iteration is given by, 



(B5) 



We have already calculated the covariance due to the first term in the previous section. We 
have to evaluate the second term and the third term. The second term is given by 
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term! = 2KL((ie {0) ) - {ei))L T K T 



(B6) 



and the third term is given by, 



term3 = KL((e {0) e {0) ) + (H) - 2(ie {0) ))L T K T 



(B7) 



Now we will evaluate the covariance the terms (e^ )^ )}, (ii), (ie^- ') . 



rn 



(^ ) 4 0) } = 7^ 0) (i + ^ ) )7f(i + 4 0) ) 

Sn) = QmpQnq(^m^n) = QmpQnq'Yp ^(1 ~\~ Hp ^)Tq (1 "I - ^) 
(^m e n') = Qmp"fp (}- + K p^)lfn "I" K n^) 



(B8) 
(B9) 
(BIO) 



Using these expressions for term2 and term3 we can evaluate the covariance in the reconstructed 
mass after two iterations. We can do a similar calculation for the next steps iteratively and keep 
higher order terms in the expansion of Equation IB2I and obtain an expression for the covariance of 
the reconstructed mass. 
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Table 2. Measuring ellipticity of dark matter and light distribution. 
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